Vortex lattice studies in CeCoIn 5 with H _L c 
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We present small angle neutron scattering studies of the vortex lattice (VL) in CeCoIns with 
magnetic fields applied parallel (H) to the antinodal [100] and nodal [110] directions. For H || [100], 
a single VL orientation is observed, while a 90° reorientation transition is found for H | [110]. For 
both field orientations and VL configurations we find a distorted hexagonal VL with an anisotropy, 
r = 2.0 ± 0.05. The VL form factor shows strong Pauli paramagnetic effects similar to what have 
previously been reported for H | [001]. At high fields, above which the upper critical field (H C 2) 
becomes a first-order transition, an increased disordering of the VL is observed. 

PACS numbers: 74.25.Op,74.25.Uv,74.70.Tx,61.05.fg 



The heavy fermion d-wave superconductor CeCoIn.5 
continues to attract great interest due to the fascinating 
interplay between unconventional superconductivity and 
magnetism exhibited by this material [IHH . In particular 
CeCoIns satisfies the stringent conditions believed to be 
necessary for the existence of a spatially modulated su- 
perconducting order parameter first described by Fuldc, 
Ferrell, Larkin and Ovchinnikov (FFLO) 0, E}: a first 
order superconducting transition at low temperatures in- 
dicating Pauli limiting and a mean free path much 
greater than the superconducting coherence length §. 
While a number of bulk measurements support a FFLO 
phase in CeCoIn 5 [H, 11 1, an experimental verifica- 



tion by a direct probe of the predicted spatial variation 
of superconductivity has not yet been reported. It is 
important because its unusual order parameter has rele- 
vance to a wide range of different areas of physics rang- 
ing from magnetars/neutron stars [b2|, u ltracold atomic 
gases [l3| to organic superconductors [IH. Furthermore, 
a field induced incommensurate antiferromagnctic order 
(Q-phasc) is observed in the same high field - low temper- 
ature region of the phase diagram where the FFLO state 
is proposed to exist [l5-17|. At present the detailed con- 
nection between the FFLO and the Q-phase remains an 
area of active research Til- 22\. 



Information about the nature of the superconducting 
state in CeCoIns can be obtained from studies of the 
vortex lattice (VL). Previous VL studies with H || [001] 
using small angle neutron scattering (SANS) have shown 
a rich vortex phase diagram and strong Pauli paramag- 
netic effects [231-427 1 . which extend well below the field 
H a where H c2 becomes a second order transition Q. 
In this field orientation, the supercurrents and quasi- 
particles travelling in the plane perpendicular to H are 



strongly affected by the d x 2_ y 2 nodes of the supercon- 
ducting order parameter and a square VL aligned with 
the nodes is formed over a significant portion of the VL 
phase diagram 25J, 127J . For H _L c, the symmetry in the 
plane _L H is no longer tetragonal and nodal effects are 
less strong, though Pauli-limiting effects are similar in 
magnitude for all field directions. 

In this letter we report the results of SANS studies of 
the VL in CeCoIns with the magnetic field applied in 
the basal plane which are the first such studies in any d- 
wave superconductor. Experiments were performed for 
field applied both along the antinodal ([100]) and nodal 
([110]) directions @, [H S-Eij]. For both field orien- 
tations a distorted hexagonal VL is observed, reflecting 
the anisotropy in the ac-plane. A single VL orientation 
is observed at all fields for H | [100], while the VL for 
H || [110] undergoes a 90° reorientation transition in the 
range 7.5 — 8.7 T. Measurement of the VL form factor 
show a strong Pauli paramagnetic effect similar to what 
we have previously reported for H || [001]. A broaden- 
ing of the VL rocking curves is observed at fields above 
Hq, suggesting that the abrupt vortex nuclcation asso- 
ciated with a first-order H C 2 leads to a freezing in of 
disorderin g. Further experimental details are provided 
elsewhere 27|, |41 1 . 



We first discuss the symmetry and orientation of the 
VL, beginning with H \\ [100]. Fig. QJa) shows a VL 
diffraction pattern with the peak positions fitted to an 
ellipse, indicating a distorted hexagonal VL. The distor- 
tion can be quantified by the VL opening angle ft. As 
shown in Fig. |Tfb) the VL structure remains unchanged 
for all measured fields and temperatures, with a constant 
opening angle, /? = 41.0°±0.7° at high fields, and perhaps 
a slight fall at low fields. The opening angle is directly 
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FIG. 1. (Color online) VL structure with H || [100] (anti- 
nodal direction), (a) Diffraction pattern obtained at 9.0 T 
and 50 mK showing a distorted hexagonal VL. Only Bragg 
peaks on the upper half of the detector were rocked through 
the Ewald sphere. The open circles show positions of Bragg 
peaks obtained by a reflection through q = 0. (b) Field de- 
pendence of the VL opening angle /? defined in (a), at 50 and 
350 mK, following either a field cool (FC) or field ramp (FR) 
procedure [4lf] . The dotted line represents H C 2 at 50 mK 



related to axial ratio of the ellipse in Fig. [TJa), which is 
a measure of the anisotropy in the plane perpendicular 
to the applied held, r 100 = tan 60°/ tan p = 1.99 ± 0.05. 
This quantity mainly reflects the anisotropy of the Fermi 
velocity in the ac-plane but includes gap effects, espe- 
cially at higher fields flisj . 

Rotating the sample such that H || [110] significantly 
changes the VL phase diagram. At low field a dis- 
torted hexagonal VL is observed, similar to that found 
for H || [100], with Bragg peaks on the major axis of 
the ellipse as shown in Fig. [2ja). However as the field 
is increased the VL undergoes a 90° reorientation tran- 
sition to a configuration with Bragg peaks on the mi- 
nor ellipse axis as shown in Fig. |2Jb). We denote these 
two VL orientations as respectively the low-field (LF) 
and high-field (HF) phases, with the former observed 
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FIG. 2. (Color online) VL with H || [110] (nodal direction). 
Diffraction pattern obtained at 90 mK and at an applied field 
of 7.5 T (a) and 9.5 T (b) respectively, (c) Field dependence 
of the VL opening angles j3, 4> and a at 90 mK. The yellow 
area indicates the transition region between the two VL ori- 
entations. The inset shows a diffraction pattern as well as a 
schematic of the rotated two-domain HF-phase. (d) VL phase 
diagram showing the extent of the LF and HF VL phases. Ho 
and H C 2 are adapted from ref. @ and Q-phase from ref. [l5|. 



for H < 7.5 T and the latter above 8.7 T. In between 
these two fields a VL rotation away from the HF-phase 
is observed as shown in the inset to Fig. [He), but is 
difficult to study systematically due to weak scattering 
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in the transition region. Fig. [2jc) shows the field de- 
pendence of the VL opening angles. The opening angle 
a, defined in the inset, remains close to 4>, suggesting 
that the transition to the LF-phase is first-order. Within 
the LF and HF phases the opening angle again remains 
essentially field-independent with (3 = 41.9° ± 1.0° and 
<t> = 31.6° ±0.5°, corresponding to r^ F = 1.94±0.07 and 
Y% = tan307tan(</>/2) = 2.05 ±0.03. While ri 00 ,r^ F 
and could be considered identical within the ex- 

perimental error, we note that the relative values are in 
qualitative a gree ment with the theoretical prediction of 
Hiragi et al. ]V^: T 100 = T\f < Tfl (although this was 
for a field of 0.211iJ c2 rj 2.4 T). 

The VL phase diagram with H \\ [110] is summa- 
rized in Fig. I2fd). Similar to what was found for 
H || [001] [25|, |27j, the transition does not coincide with 
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the appearance of the magnetic Q-phase [15H17} nor the 
proposed FFLO-phase ((Ha [13, El • More importantly, 
our results demonstrate a distinct difference between the 
antinodal and nodal directions in CeCoIns. In contrast to 
the VL anisotropy, T, our phase diagram does not agree 
with the theoretical predictions of opposite VL orien- 
tation transitions with increasing temperature for fields 
along [100] and [110] respectively 18|- This is not sur- 
prising since calculations of the VL structure depend sen- 
sitively on the Fermi surface and gap anisotropy and are 
notoriously difficult due to the small energy differences 
involved. This illustrates the need to include realistic 
material parameters in such calculations. 

It is interesting to compare the VL behavior to that 
reported for LuNi 2 B 2 C and ErNi 2 B 2 C with H || [100], 
where a 90° VL first order reorientation is observed [3l| . 
Thermal conductivity measurements in non-magnetic 
LuNi 2 B 2 C and YNi 2 B 2 C have suggested the presence of 
point-nodes along the (100) directions 32j,l33(, indicating 
that the VL reorientation is observed in cases where the 
field is applied along the nodal direction and is unrelated 
to a vortex core expansion which is believed responsible 
for the reentrance of the square VL phase in CcCoIns 
with H || [001] H3. 

We now turn to studies of the VL form factor F(q), 
which is a measure of the magnetic field modulation in 
the mixed state. To avoid complications arising from a 
changing VL orientation, and the difficulty obtaining a 
well-ordered single domain VL in the transition region, 
these measurements were performed only for H || [100]. 
The form factor is related to the integrated reflectivity R 
of the VL, which is obtained by rotating and/or tilting 
the cryomagnet and sample such that the VL scatter- 
ing vectors cut through the Ewald sphere. Examples of 
rocking curves obtained in this fashion are shown in the 
inset to Fig. @] Gaussian fits to the rocking curves give 
an integrated intensity which is divided by the incident 
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FIG. 3. (Color online) Field dependence of VL form factor, 
(a) Form factor with H || [100]. The dashed lines are guides 
to the eye and Hdi is shown for T — 50 mK. (b) Comparison 
of |F(q)| 2 measured with field along [100] and [001] [27| . (c) 
Multiplying the form factor for H \\ [100] by (2.55) 2 = 6.5 
causes |F(g)| 2 for the two field orientations to fall on a single 
curve. 



neutron flux to yield the integrated reflectivity 



R = 



in 



where 7 = 1.913 is the magnetic moment of the neutron 
in nuclear magnetons, A n is the neutron wavelength, t is 
the sample thickness, </>o = h/2e — 2067 T nm 2 is the 
flux quantum and q is the length of the VL scattering 
vector. The resulting VL form factor for H \\ [100] is 
shown in Fig. [3ta). Contrary to the commonly observed 
exponential decrease with increasing field [34[ , |-F(<7)| 2 in- 
creases with increasing field up to a maximum at 10 T 
after which it decreases on approaching the upper critical 
field. We note that a qualitatively similar behavior is ob- 
served at 50 and 350 mK and for different VL preparation 
(FC/FR). However there are subtle differences between 
the form factor above Hq = 10 T when one compares 
VLs obtained following a FC and a FR. We will return 
to this point later. 

The unusual field dependence of the form factor is sim- 
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ilar to what has been reported for TmNi2B2C [35j and 
CeCoIn 5 with H || [001] HH^, and can at least 
partly be explained by Pauli paramagnetic effects [35l — 
|37J. Briefly, there is a significant spin-polarization of 
the unpaired quasi-particles in the vortex cores, result- 
ing in an increased amplitude of the modulation of the 
magnetic field. Moreover, it was recently shown that ac- 
counting for antiferromagnetic fluctuation in the vortex 
cores yields both a skewing and an increase of the form 
factor |39j|, in better agreement with our findings. 
One notes that while the maximum in \F(q)\ 2 coincides 
with the onset of the Q-phase at 50 mK, the two effects 
do not appear to be linked. This is supported by the 
weak T-dependence of the field at which the form fac- 
tor is maximum and that a similar behavior is observed 
for H || [001] [27J where no Q-phase has been observed. 
Moreover, the temperature dependence of the VL peak 
intensity (top of rocking curve) at 10.9 T showed no fea- 
tures except an abrupt drop to zero at H C 2- 

Fig. (3Jb) and (c) shows a comparison of the form fac- 
tors obtained with the applied field perpendicular and 
parallel to the crystalline basal plane as a function of 
the reduced field [40J , which can be collapsed onto a sin- 
gle curve by multiplying Fioo by a constant factor of 
2.55 ± 0.1. As we have previously shown j^, the form 
factor at H C 2 is directly related to the jump in the mag- 
netization, yielding AM001/AM100 ~ 2.8 [3[ which is in 
quite good agreement with the scaling of |-Fioo|- What 
is remarkable is that the rescaled form factors agree as a 
function of reduced field over the entire measured range, 
whereas one would expect the ratio to gradually approach 
the Fermi velocity anisotropy (~ 2) at low fields where 
Pauli paramagnetic effects are small. This scaling invites 
detailed theoretical analysis. 

Finally we return to the observation that above Hq 
the VL acquires a field/temperature history dependence, 
leading to a difference between the FC or a FR case. As 
shown in Fig.[3ta), the measured form factor is lower fol- 
lowing a FC than following a FR for fields above 10 T, 
indicating a more well-formed VL in the latter case. The 
same trend is seen in Fig. @] showing the field depen- 
dence of the rocking curve widths, which is inversely 
proportional to the longitudinal correlation length of the 
VL. Below 10 T the width is resolution limited. Above 
Hq, one sees a striking increase in the FC rocking curve 
widths. We speculate that this is due to the abrupt 
appearance of the vortices while slowly cooling through 
the first-order superconducting transition, immediately 
freezing in of any disordering. Furthermore, the VL can 
be effectively annealed by inducing vortex motion (FR), 
as shown in the inset in Fig. @] which compares FC and 
FR rocking curves for 10.5 T. 

In summary we have studied the VL in CeCoIns with 
H _L c. While a single, distorted hexagonal VL is found 
for H || [100], a 90° reorientation is observed for H \ 
[110]. The VL form factor shows Pauli paramagnetic 
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FIG. 4. (Color online) VL rocking curve widths correspond- 
ing to the data in Fig 1(b). The significant scatter in the 
data at low fields is due to the low scattered intensity and 
higher background at shorter q. The calculated experimental 
resolution is shown by the dashed line. The insert compares 
the rocking curves obtained at 10.5 T following a FC and FR. 
The lines are Gaussian fits to the data and 60 denotes the 
center of the rocking curve. 



effects, similar to what is observed for H || [001]. Finally 
both the form factor and the rocking curve widths show 
a clear dependence on the VL preparation above Hq. 
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